SU-ITP-93-29 
IEM-FT-77/93 
|astro-ph/9312039| 
16 December 1993 



Fluctuations of the Gravitational Constant 
in the Inflationary Brans— Dicke Cosmology 



O ■ Juan Garcia— Bellido 1 , 

O 

3 



Andrei Linde 2 , 

Department of Physics, Stanford University, Stanford, CA 94305-4060, USA 



^ ! Dmitri Linde 3 

OV 

CO . California Institute of Technology, Pasadena, CA 91125, USA 

O 
<N 

CO 

Abstract 

According to the Brans-Dicke theory, the value of the gravitational constant G which we 
measure at present is determined by the value of the Brans-Dicke scalar field at the end of 
inflation. However, due to quantum fluctuations of the scalar fields produced during inflation, 
the gravitational constant G(4>) may take different values in different exponentially large parts of 
the Universe. We investigate the probability distribution P p to find a domain of a given volume 
with a given value of the gravitational constant Gat a given time. The investigation is performed 
for a wide class of effective potentials of the scalar field a which drives inflation, and with two 
different time parametrizations. Our work is based on the analytical study of the diffusion 
equations for P p , as well as on the computer simulation of stochastic processes in the inflationary 
Universe. We have found that in some inflationary models the probability distribution P p rapidly 
approaches a stationary regime. The shape of the distribution depends, however, on the choice 
of the time parametrization. In some other models the distribution P p is not stationary. An 
interpretation of our results and of all ambiguities involved is outlined, and a possible role of 
anthropic considerations in determination of the gravitational constant is discussed. 
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1 Introduction 



One of the most amazing properties of inflationary cosmology is the process of self-reproduction 
of inflationary domains of the Universe (for a review see PQ). This process exists in many versions 
of inflationary universe scenario, including old inflation [2] , new inflation [3J 0] , chaotic inflation 
[5] and extended inflation Self-reproduction of inflationary domains implies that there is no 
end of the evolution of the Universe. This process divides the Universe into many exponentially 
large domains with all types of symmetry breaking and with all types of compactification of 
space-time compatible with inflation. 

The best way to describe the global structure of the Universe in this scenario is provided 
by the stochastic approach to inflation. The original version of this approach [7j was based 
on the investigation of the distribution of probability P c (<f),t) to find a given field at a given 
time at a given point. This approach was not well suited for the investigation of the process of 
self-reproduction of the Universe. A more adequate approach is based on investigation of the 
distribution of probability P p ((j),t) to find a given field at a given time in a given physical 
volume El El- 

The most detailed study of the distribution P p was performed recently in ^H] , where it was 
shown that in many inflationary models including the models with the effective potentials <p n and 
e°"^ this probability distribution rapidly approaches a stationary regime. This means that if one 
takes a section of the Universe at a given time t and calculates the relative fraction of domains of 
the Universe with given properties (with given density, with given values of various scalar fields, 
etc.), the result will not depend on the time t, both during inflation and after it. 

This result represents a major deviation of inflationary cosmology from the standard Big 
Bang paradigm. A lot of work is still needed to verify this result and to obtain its consistent 
interpretation in the context of quantum cosmology. One should also study how our methods work 
in the context of more complicated models, including several different scalar fields [TT) IT21 ITUj. 

A natural idea would be to apply our methods to extended inflation scenario |T3|, which is 
a hybrid of the Brans-Dicke theory and old inflation. In this scenario one has two scalar fields: 
the Brans-Dicke field <fi and the inflaton field a. However, during extended inflation only the 
former evolves in time, which reduces the problem to the one we have studied already. Moreover, 
the value of the scalar field (f) after inflation in this theory typically is determined not by the 
stochastic dynamics during inflation, but either by the position of a pole of a non-minimal kinetic 
term of the field <fi or by the effective potential of this field which should be introduced into the 
theory in order to make it consistent with observational data [14] . 

In this paper we will consider theories of another type, which are hybrids of the Brans-Dicke 
theory with new ^SlCSl or chaotic inflation (TJj. In these models, especially in the chaotic inflation 
model of ref. [T2j, simultaneous stochastic evolution of the fields </> and a is very nontrivial. One 
of the most interesting consequences of this evolution is the division of the Universe after inflation 
into exponentially large regions with different values of the gravitational constant in each of them. 
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Indeed, the gravitational constant G in the Brans-Dicke theory is a function of the field 0, 

c =«^- (i) 



where u is the Brans-Dicke parameter, u > 500 [T8j[T9]. The value of the scalar field practically 
does not change after inflation. Therefore, investigation of the distribution P p (<p,t) gives us the 
fraction of the physical volume of the Universe where the gravitational constant G has the effective 
value ((TJ). 

The stationary character of the probability distribution P p (<j) in the theory of one scalar field a 
in the standard general relativity theory is closely related to the existence of the Planck boundary 
cr p , where the potential energy density V(o~) becomes comparable with the Planck density M£. 
Typically the distribution P p (a,t) rapidly moves towards large a, for the reason that the volume 
of domains with large V(o~) grows very fast. The distribution P p (a,t) becomes stabilized as it 
approaches the Planck boundary, where, as it is argued in [TOj, the process of self-reproduction 
of inflationary domains is less efficient. 

In the Brans-Dicke theory the situation is somewhat different. The Planck boundary is not 
a point, but a line <p p (a p ), where 

4vr 2 

V(op) = M p 4 (0 p ) = — 0* . (2) 

Therefore, after the distribution P p (a, 0; t) approaches the Planck boundary, it can still move 
along this boundary. One may encounter three different possibilities: 

1) The distribution P p (a, 0; t) never becomes stationary for any values of and a. Nevertheless 
the properties of domains of inflationary universe filled by the fields and o do not depend on the 
time when these domains were formed. This is the most profound stationarity which is present in 
all inflationary models where the process of self-reproduction of inflationary domains is possible. 
In ref. [TU] we called this property micro stationarity, or local stationarity. 

2) The distribution P p (a,4>;t) normalized over all possible values of and a gradually ap- 
proaches a stationary regime. In ref. ^U] we called this property macro stationarity, or global 
stationarity. 

3) The maximum of the distribution P p (a, 0; t) runs away to infinitely large (or to infinitely 
small) values of and a. As a result, the relative fraction of volume filled by any finite and 
cr as compared with the total volume of the Universe decreases in time, which means that the 
global stationarity is absent. However, one may wish to exclude from consideration domains with 
infinitely large (or infinitely small) values of and a and concentrate on some finite part of space 
(0, cr) (for example, on those values of and a which are consistent with the existence of life as 
we know it). Then in some cases one may find out that the probability distribution P p (a,4>;t) 
normalized over this part of space (0, a) (or the ratio prf~|~^ ) approaches a stationary regime. 
In this case we will speak of a runaway stationarity. 

As we will see, all these possibilities may be realized in the inflationary Brans-Dicke cosmology, 
depending on the choice of the effective potential V(a). 
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There exists another important issue to be discussed in this paper. Even if the distribution 
P p (a, 0; t) does not depend on time, it may depend on the choice of possible time parametrizations 
[TUj . As we will see, in the Brans-Dicke theory this dependence can be quite strong. 

In Section El of this paper we describe chaotic inflation in the Brans-Dicke theory at the 
classical level, following ref. [T7] . In Section |3] we derive the stochastic equations describing 
the evolution of quantum fluctuations of the scalar fields and a during inflation, as well as 
the origin of adiabatic energy density perturbations from these fluctuations. In Section 0] we 
describe the phenomenon of self-reproduction of the inflationary universe in the presence of 
both scalar fields as well as the stochastic approach to inflation in Brans-Dicke theory. We 
study the boundary conditions for P p in Section These equations are extremely complicated, 
and it is not always possible to solve them analytically. Therefore we performed a computer 
simulation of the stochastic evolution of the scalar fields during inflation. We describe these 
numerical simulations and their results in Section H In Section [7| we analyze the case of a 
constant vacuum energy density, where runaway stationarity appears. For general increasing 
potentials V(a), the probability distribution in the physical frame will rapidly move to large 
values of a, constituting what we called runaway solutions. These non-stationary solutions are 
described in Section |H| Imposing boundary conditions at large o or introducing steep potentials 
V(a) will give a stationary distribution, as described both numerically an analytically in Section 
IH1 In Section El we derive the stochastic equations using a different time parametrization, where 
instead of the usual time t we choose the time r ~ lna(t) [TU]. The methods of computer 
simulations of stochastic evolution in time r are different from the methods which we use in our 
simulations of the evolution in time t. We describe these methods and their results in Section 
ITT1 We compare the results obtained in different time parametrizations and make an attempt of 
their interpretation in Section 

As we already mentioned, under certain conditions the probability distribution P p is station- 
ary, i.e. time-independent. This apparently gives us a possibility to calculate the position of the 
maximum of the distribution P p and thus predict the most probable value of the gravitational 
constant in the Universe. However, this idea involves some ambiguous speculations. The readers 
will find them in the Appendix. 



2 Inflationary Brans-Dicke Cosmology 



In this section we will describe the classical evolution of the inflaton field with a generic chaotic 
potential, in the context of the Jordan-Brans-Dicke theory of gravity, 



S = J d^Xy/^g 



107T Z Z 



(3) 



Here M 2 (<f>) = ^0 2 , <f) 2 /8uj is the Brans-Dicke field, and a is the inflaton field. We will consider 
several different potentials, including V(a) ~ ^a 2 , \a A , -^(m 2 — \a 2 ) 2 and ^a 2 + |cx 4 log^. 
The theory (JHJ) looks similar to the extended inflation model ^3]- The difference, which will be 
very important for us, is that in the extended inflation scenario inflation occurs at a — 0, whereas 
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in our case inflation occurs during the slow rolling of the field a 17 . The complete equations of 
motion in a FRW Universe with scalar fields <fi and a are 

D 2 + ±.r) <f>(t) = , 

D*a(t) = -q; , (4) 



2 « 4a; /l - 2 1 . 2 



where = ±1,0; £> 2 = d 2 /dt 2 + 3Hd/dt — k 2 /a 2 , and if = a/a. During inflation we can 
consistently use the slow roll-over approximations, V"(a), k 2 /a 2 -C H 2 , <fi <C Hcf) <C H 2 <p, 
\a 2 + \(f) 2 <C V(p). The equations of motion (jlj) then simplify to 

V'(a) 



-H 



a 



uj 3H ^ 

R = -12H 2 , H 2 = ^V{a) . 

We will be most interested in theories with V(a) = ^cr 2n , for which equations (jSj) reduce to 



a \6u J 



It follows that in those theories <\> and cp = y ^ a move along a circle of constant radius in 
the plane (ip,4>). We can parametrize the classical trajectory by polar coordinates (0(t), f{t)) = 
(r sin 8 (t),r cos 8 (t)) with constant r, and angular velocity 

8\t)=( f] 1/2 (^[r cos 8(t)r^ (7) 



3u) \2 

For n = 1 we find solutions with a constant angular velocity 

(pit) = r sin O + -=t . 



r / m 



= cos \0 o + -j=t) , (g) 



sin 6 1 , 



a(i) = a Q 



Here < 8 Q < tt/2. These solutions in the interval 8 Q : ^ L < t ^ corres P on d to the usual 
power-law behavior a(t) ~ t u . For n > 2, the angular velocity decreases with time and the 
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classical solutions are more complicated, but inflation is still power-law. These solutions are 
actually attractors of the complete equations of motion (jljl. for all n |2(Jj . 

We will now discuss initial and final conditions for the inflationary universe. In the chaotic 
inflation scenario, the most natural initial conditions for inflation are set at the Planck boundary, 
V(cr p ) = Mp(0 p ), beyond which a classical space-time has no meaning and the energy gradient of 
the inhomogeneities produced during inflation becomes greater than the potential energy density, 
thus preventing inflation itself. The initial conditions for inflation are thus defined at the curve 

/ \ 1/2 / . \ 1/4 , 

( t ) p = (|^) (2^) a p ■ On ^ ne other hand, inflation will end when the kinetic energy density 
of the scalar fields becomes comparable with the potential energy density, |0 2 + \a 2 ~ V(a) 
or m 2 = V"(a) = H 2 ((p,cr). This condition corresponds to the end of inflation boundary <j e = 



We will also consider an inflaton potential which leads to spontaneous symmetry breaking, 
V(a) = | (cr 2 — cr 2 ) 2 , where cr D = In this theory equations (jSJ) are expressed as 



a 



(9) 




The fields move approximately along a circle centered at a = a Q , in clockwise direction if the 
initial condition is to the left of the minimum of the potential and in anticlockwise direction if it 
is to the right. In fact, for a ^> a they will soon approach the asymptotic solution to © with 
n = 2, while for a <^ a , we find the solution 
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m 



a(t) 
a(t) 




(10) 



Inflation will now occur in two disjoint sectors, either to the left or to the right of the minimum. 
The Planck boundary for initial conditions is again defined by V(cr) = Mp (</>), while the end of 
inflation is given by the condition |cr 2 ~ V(a). These two boundaries are defined by the curves 




°o ~ (J p\' , 



In the absence of any potential for <p, the Brans-Dicke field remains almost constant after 
inflation, and therefore the Planck mass today is approximately given by its value at the end 
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of inflation, M p ~ J ^ <p e . On the other hand, the total amount of inflation is approximately 

a p \<t>p) 

The value of the Brans-Dicke parameter u is bounded by post-Newtonian experiments [T%] 
and primordial nucleosynthesis [19J to be very large, u > 500, and therefore it is appropriate to 
use the approximation u 3> 1 in the following analysis. 



3 Stochastic Methods 



In this section we will describe the stochastic evolution of the scalar fields during inflation and 
the generation of adiabatic energy density fluctuations from quantum fluctuations of these fields. 



3.1 Quantum Fluctuations of Scalar Fields 

The classical scalar fields (j) and a in de Sitter space are perturbed by their own quantum fluctu- 
ations, which are stretched beyond the horizon and act on the quasi-homogeneous background 
fields like a stochastic force. In order to calculate this effect we will coarse-grain over a horizon 
distance and split the scalar fields into long-wavelength classical background fields (p(x) and a(x) 
plus short-wavelength quantum fluctuations with physical momenta k/a > H, 



(12) 



where e is an arbitrary parameter that shifts the scale for coarse-graining |7 . The physical results 
turn out to be independent of the choice of e. The quantum fluctuations are assumed to satisfy 
the following commutation relations 



4>(x,t) = 4>(x,t) + j d 3 k 9{k — eaH) a k Uk(x) + a\u* k {x) 
a(x,t) = a{x,t) + / d 3 k' 9(k' - eaH) \byVy{x) + b\,v* k ,(x) 



[a k , al] = [h, b\,] = 5 3 (k - k') , [a k , b\,} = . 



(13) 



The exact solutions to the scalar fields' equations (jU) in de Sitter space with V(a) = \m 2 a 2 
are given by |2*Tj 



gifc-x Hti 9 3 9 

«*(*)= J^—V^H^(k V ), ^=- + -^- ) 

/ \ H V i — iT(2)n> \ 2 9 m 2 9 



( 2vr )3/2 2 v ' " v " ' 4 H 2 4 ' 

where r\ = — (aH)' 1 is the conformal time, and H^ 2 (x) = — ^/^exp(zx) (l + ^j. The amplitude 



S 



of the quantum fluctuations of and a can be computed as 



i i-U |2\V2 H 



5a = [A'irk 6 \v k 



/3| |2\V2 # 



27r (15) 



2vr ' 

which coincides with the Gibbons-Hawking temperature. 

These quantum fluctuations then act as a stochastic force on the classical background fields. 
One could write the evolution of the coarse-grained fields in the form of Langevin equations 

8(f) M 2 {<j))dH , H 3 ' 2 



dt 2n d(f) 2n 

da MjffldH H*'\ 1 ) 

dt " An da + In U } ' 

where £ and £ behave like an effective white noise generated by quantum fluctuations, (((t) £(£')) 
= = — an d (C(0 £(0) = 0; which leads to a Brownian motion of the classical 

scalar fields and a, with a typical step (fT5|). 

We can alternatively describe the stochastic process in terms of the probability distribution 
P c (a, 0; t). This distribution describes the probability to find the fields and a at a given time t 
in a given point. Equivalently, it describes the probability to find, at a given time t in the domain 
with a given comoving (i.e. non-expanding) volume, the fields with mean values and a. As it 
was shown by Starobinsky this probability distribution satisfies the Fokker-Planck equation, 

dP c _ 8 (Ml{<P)dH H 3 / 2 d{H 3 / 2 P c )\ 

' Pc 




dt da \ An da 8n 2 da J , , 

Mj(<j>) dH H 3 l 2 d{H 3 ' 2 P c ) \ = _dJ^ _ d^ [L7) 
" 1 2n 90 8n 2 90 J da 8(f) ' 

where we have chosen the Stratonovich version of stochastic processes. This equation can be in- 

dP c 

terpreted as the continuity equation — — h V • J = associated with the conservation of prob- 

dt 

ability. The first terms of each current correspond to the classical drift forces for the fields 
and a (jBJ), while the second terms correspond to the quantum diffusion due to short-wavelength 
fluctuations (fTHJl . 

One can then compute the field correlations during inflation with the help of this probability 
distribution, assuming that H is approximately constant, 

| <0V> = ^- 2 p( P - l) (0»-V) + i| 

+ ^q(q ~ 1) - X (0V- 
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For example, for V(a) = \m o , we find 

(**>= ^(^-o-^. 

2m 2 1 \ 

#3 #2 v-; 

Note that the dispersion of both fields due to Brownian motion is identical for a relatively large 
time interval. 

We will now describe the origin of energy density perturbations in the Universe from quantum 
fluctuations of the scalar fields during inflation, and postpone the study of the self-reproduction 
of the inflationary universe for the next section. 



3.2 Energy Density Perturbations 

In this subsection we will describe the generation of adiabatic energy density perturbations, 
when both fields a and <f> are included. The amplitude of these perturbations has been computed 
previously in [221 120] • We will give here an alternative derivation. 

The energy density perturbations could have originated during inflation in our model ©, as 
quantum fluctuations of both scalar fields that first left the horizon during inflation and later 
reentered during the radiation or matter dominated eras. The amplitude of those perturbations 
can be computed in the Einstein frame {g^ u = 4> 2 9^u) by using the equality [23J 



p \ 8p 



3p + pl p 



2HC 



3p+pl p 



(20) 

IHC 



where IHC corresponds to the time when the perturbations first left the horizon and 2HC to 
the time when they reentered. Since during inflation p + p — 0, the amplitude of reentering 
perturbations can be written as 



5_p 
P 



2rj Sp 

2HC~ 3 P+P 



(21) 

IHC 



where rj = 2/3 (3/5) for perturbations reentering during the radiation (matter) era. During infla- 
tion, the pressure and energy density in our theory (j3J) in the Jordan frame have the expressions 

P = U 2 + \° 2 + V(a) , 

2 2 (22) 

P = + - V{a) , 

and thus p + p = <fi 2 + a 2 . The adiabatic energy density perturbations follow from the quantum 
fluctuations of the fields. In the Einstein frame, y/gp = y/gp, 

p V{&) 
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Using the equations of motion (JHJ) we find (for the cold matter dominated Universe) 



Sp 
P 



2HC 



- 6 -H 

5 



<j)5<j) + &Sa 



+ a 



V2 



3H 2 



a 



IHC 



57T 02 _|_ o-2 



(24) 



where N\ stands for the number of e-folds before the end of inflation associated with the horizon 
crossing of a particular wavelength. For perturbations of the size of the present horizon, we must 
compute the last expression for N\ ~ 65 1 . For a general potential V(a), we can write the 
density perturbations as 



Sp 
P 



1 + 



24#(<x,0) V{a) 
T M p 2(0) V'{a) VT + W 



(25) 



Note that in the large u limit, cf) <C a during the last stages of inflation, which ensures the 
approximate equivalence of the Einstein and Jordan frames. We then recover the usual expression 

m 



Sp 

p 



24H(a,(p) V(a) 



5 Af2(0) V'{a) 



(26) 



where M p is now 



-dependent. For theories with potentials of the type Xa , it behaves like 



Sp 
P 



6u ( 2ujX 



V2 a n+l 



5mr \ 3n 



(27) 



Af=65 



In the case of the theory Act 4 , the density perturbation (|2T|) takes the usual vA dependence. 
However, for the theory m 2 a 2 /2, the perturbation on the horizon scale is given by 



Sp 
P 



50 m 

MJd>: 



(28) 



Therefore, we note that the larger is the Planck mass at the end of inflation in a given region of 
the Universe, the smaller will be the density perturbation in this region for this model. We will 
return to the discussion of this result at the end of the paper. 



4 Self— reproduction of the Inflationary Universe 

Quantum perturbations produced during inflation are responsible not only for galaxy formation, 
but also for the process of self-reproduction of the whole inflationary universe. This is the 
important effect which we are going to consider here in the context of Brans-Dicke cosmology. 
We will begin with an elementary description of this effect, and then return to its description in 
the context of the stochastic approach to inflation. 
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4.1 Elementary Considerations 



An important property of the inflationary universe is that processes separated by distances I > 
H^ 1 proceed independently of one another pQ. In this sense any inflationary domain of initial 
radius exceeding if -1 can be considered as a separate mini-universe, expanding independently of 
what occurs outside it, as a consequence of the "no-hair" theorem for de Sitter space. 

During a typical time if -1 each such domain expands e times, and its volume grows e 3 ~ 20 
times. This means that this domain becomes divided into 20 independent inflationary domains. 
The values of the classical fields <fi and a inside each of these domains can be obtained by solving 
classical equations of motion for these fields. However, in addition to it one should take into 
account quantum fluctuations of these fields, which become "frozen" during the time H~ x . They 
have a typical amplitude but they may have different signs in each of the new 20 domains. 
Let us assume that this amplitude is much greater than the classical shift of the values of these 
fields during the time H~ x . In this case in 5 out of 20 domains the field jumps towards its 
smaller values, and the field a jumps towards its greater values. The same happens during the 
next time H~ l . 

In the context of theories with V(a) ~ Acr 2n , this leads to a continuous process of recreation 
of inflationary domains. Using the classical equations of motion (jUJ), one finds that the condition 
that quantum diffusion is more important than the classical drift of the fields <fi and a is satisfied 



where the last term corresponds to the Planck boundary. From this equation it follows that the 
inflationary universe with most natural initial conditions (i.e. not far from the Planck boundary) 
enters eternal regime of self-reproduction. During this regime the Universe becomes filled with 
all possible values of the fields <fi and a, independently of their initial values in the region (J29|) . 
Consequently, the parts of the Universe where inflation ends will consist of many exponentially 
large domains in which the Planck mass M p (0) and the gravitational constant G(4>) may take all 
possible values from to oo. 

This scenario, which we called eternal inflation j3], deviates strongly from the standard Big 
Bang theory. For example, according to the standard theory a closed Universe should eventually 
collapse and disappear. In our scenario a closed Universe which has at least one inflationary 
domain with a field <fi in the interval (j2Uj) will never disappear as a whole. 

Here one should make some comments to avoid terminological misunderstandings which some- 
times appear in the literature. Eternal inflation does not mean that each part of the Universe 
eternally inflates. The typical length of each geodesic at the stage of inflation is finite. However, 
there is no upper limit to the length of these geodesies, and those extremely rare geodesies which 
have large length give the dominant (and permanently growing) contribution to the total volume 
of the Universe. Therefore in our scenario inflation in the whole Universe has no end, even though 
it ends on each particular geodesic within a finite time. 

One may try to reverse the question and ask whether inflationary universe has any beginning. 



for UHEi 




(29) 
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Unfortunately, the answer to this question is much less definite. One may argue that each geodesic 
being continued to the past has finite length (it begins with a singularity) j2H|. However, this 
is not enough to prove that the Universe has a single beginning at some moment t = in the 
past (Big Bang). Whereas such a possibility is not excluded, in order to prove it one should 
show that there is an upper limit to the length of all geodesies continued to the past. Indeed, 
even if long geodesies are extremely rare, they may give exponentially large contribution to the 
present volume of the Universe. At present we do not have any proof that that there exists any 
upper limit to the length of all geodesies continued to the past. Therefore finiteness of length of 
each geodesic being continued to the past does not mean yet that inflation is eternal only in 
future. We emphasize again that the length of each particular geodesic at the stage of inflation 
is also finite, and still we are speaking about eternal inflation. 

On the other hand, the properties of each particular inflationary domain created in the process 
of self- reproduction of the Universe do not depend on the time when it was created (microsta- 
tionarity). Therefore by local observations which we can make inside our domain we cannot come 
to any conclusion about the time when the Big Bang happened. Therefore, our scenario removes 
the Big Bang to the indefinite past and in this sense makes its possible existence almost irrelevant 
|lUj . In particular, the stationary probability distributions which we are going to obtain will not 
depend on initial conditions at the beginning of our computer simulations. 



4.2 Stochastic Approach 

A formal method to describe the process of self-reproduction of inflationary domains is given by 
the stochastic approach to inflation. One of the possibilities is to go along the lines of ref . jH] , to 
solve the Fokker-Planck equation (|T7j) for the distribution P c , and then to study P p using these 
solutions. 

For initial conditions of the fields a and </> far away from the Planck boundary, the proba- 
bility distribution in the comoving frame behaves like a Gaussian centered around the classical 
trajectory (cx c (t), 4> c {t)) in the (cr, </>) plane, 

f (v-vc(t)) 2 (<p-<p c (t)) 2 \ , s 



with dispersion coefficients [2 

A 2{+\ A ( U \ °c / 4 4\ A ( U \ °c °o 




3n 2 \2irJ <ft v ° CJ 3n 2 V2tt ) ^ 



(31) 



These results are very similar to the results of the investigation of P c in the theory of a single 
scalar field obtained in [Hj. One may then use the fact that during small time intervals At the 
probability distribution P p (a, (ft; At), which takes into account the difference of the rates of the 



13 



quasi-exponential growth of the proper volume in different parts of the domain, is related to P c 
in a rather simple way, 

P p (a, <f>; At) ~ P c (a, 0; At) e 3HAt . (32) 

With the help of this approximate relation one can study the qualitative features of the 
behavior of P p and confirm the existence of the regime of self-reproduction jS]. However, to 
obtain a more detailed information about P p one should study directly the diffusion equation for 
P p . This equation differs from the equation for P c only by the presence of an extra term 3HP P 

pun]: 

dP± = d (MjMdH d (H3/2p 

dt da \ 4tt da p 8tt 2 da V % - , 

Apart from studying the distribution of fields a and </> m all domains during inflation, we will 
calculate the volume of all domains where inflation ends in a state with given <fi within each new 
time interval. This gives us the fraction of the volume of the Universe where inflation ends at a 
given time t within a given interval of values of the field <fi. We call this new distribution V p (4> e , t). 
This distribution is closely related to P p . For example, in the theories with V(a) = ^a 2n 

P p (<j>e,t) ~<%-P p (<j> e ,<T e ,t) . (34) 

Indeed, during the time At all domains in the interval Act from a e — a At to a e will cross the 
boundary of the end of inflation at a = a e . According to eq. (JUJ), in the theories we consider 

JL / ,\l/2 

Aa = -a At = f f a n At. The value of the field near the end of inflation almost does not 

change, <fi = e , and a = a e = -§=0 e - This yields Act = J - (R 7 C/ 2 4>e^- Obviously, the fraction 





of the volume of the Universe where inflation ends at a given time t within a given interval of 
values of the field = e is proportional to P p (4> e , &e, t)Aa. This gives eq. (|34|). up to an overall 
normalization factor. 

Since the value of the effective Planck mass M p (0) after inflation practically does not change, 
this distribution is most directly related to the fraction of the volume of the post-inflationary 
universe with the Planck mass M p (0 e ) ~ J^- (f> e . In this paper we will be interested mainly 
in stationary distributions. Whenever our distributions will be time-independent, we will write 
them simply as P p {4>) or V p ((f> e ). 



5 Stationary Probability Distributions 

In general it is very difficult to find any analytic solutions to the equation (|3!?|) for P p (a,cj);t). 
However, in certain cases the corresponding solutions in the limit of large t can be represented in 
the simple form, 

F,(^;t)~e a P p (^), (35) 
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where E is some constant [HHHI]. In such cases the normalized probability distribution P p (a, 0) 
will be stationary. 4 Analytical investigation of P p (cr, 0) often can be simplified if one studies 
instead the function \l/(er, 0), where 

P p (a, 0) oc H- 3 ' 2 (<t, 0) exp {^^) <P) ■ (36) 

Using the identity D 2 g — 2D(Df g) = D 2 (e~f g) — [(Df) 2 + D 2 f]g, one can show that the new 
function \l/(<7, 0) satisfies a two-dimensional Schodinger-like equation 

(h 3/2 ^A 2 * + (# 3/2 ^) 2 * - ^fa = ( 37 ) 

with a new effective potential 



V(a, 0) = V, 0)V» 2 + ^ V, *) ~ V'V) 

647T 4 fi7T 2 

+ ^— FT 3 (a, 0)V(a) + —H(a, 0) - 24tt 2 #( ( 7, 0) . 

Six) UJ 



(3* 



This equation (or equation (f53J)) should be supplemented with boundary conditions. There 
are three possible boundaries in the (a, 0) plane. 

1. End of inflation boundary. Our diffusion equations are valid only during inflation. 
Therefore some boundary conditions should be imposed at the boundary where inflation ends. 
These conditions follow from the continuity of the probability distribution P c and of the prob- 
ability current J [TO]. In the theories with V(a) ~ o 2n and ui ^> 1 the field at the end of 
inflation almost does not change. The continuity condition can be expressed in terms of the field 
a changing from the right side of the boundary <j e (from a e +) to the left of it (to cr e -): 

P c (a e+ ) = P c (a e -) , J(a e +) = J(a e -) . (39) 

As it is shown in [TO], this leads to the following boundary condition on P p : 

^(if 3 / 2 (a,0)P p ( ( T,0))| cnd = O. (40) 
One can re-express this boundary condition in terms of the redefined function \1/ (cr, 0) (|36j) as 



oa \ 



= . (41) 

end 



2. Planck boundary. The distribution P p (cr, 0;t) typically tends to be shifted towards the 
region of greatest possible Hubble constant, which ensures exponentially fast growth of volume 

Since the difference between P p and P p is only in the normalization, we will usually omit the tilde and write 
P p simply as P p . 
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of inflationary domains ~ e 3H(a,<f>)t ^ However, one may argue that inflation destroys itself at 
values of the potential energy density above the Planck scale by production of large gradients of 
density. Furthermore, the classical space-time in which inflation takes place cease to make sense 
above the Planck scale, where quantum fluctuations of the metric are important. Therefore it 
is natural to impose some boundary conditions at the Planck boundary which would not allow 
a nonvanishing P p at densities higher than M p (<p). As it is argued in [TUj . most of the results 
are not very sensitive to a particular choice of such boundary conditions (absorbing, reflecting, 
etc.). Therefore we will simply assume that the probability distribution P p (a, 0; t) vanishes when 
l»=M p 4 (0): 

P p (a p , P ; t) oc ^(cr p , 4> p ) = . (42) 

Here a p , (f) p is any pair of values of the fields a and belonging to the line V(a) = M p (0) (Planck 
boundary). 

3. Boundary at large a. The two boundary conditions mentioned above are not enough 
to ensure stationarity of solutions, since the maximum of the probability distribution may move 
along the Planck boundary. The reason is very simple. In the ordinary inflationary theory the 
maximum of the probability distribution moves towards the Planck boundary since near this 
boundary the rate of exponential expansion of the Universe is maximal. In our case the Planck 
boundary is not a point where V(a) = M p but a line 0(a) (J2J). The greater is along this line, 
the greater is the energy density there, the greater is the rate of expansion. Therefore one may 
expect the probability distribution P p to move along the Planck boundary towards greater and 
greater values of and a. 

This would mean that there is no macrostationarity (global stationarity) in our model, whereas 
the microstationarity (local stationarity) is still present. Even though greater and greater num- 
bers of inflationary domains will contain indefinitely large values of the fields, there will be 
exponentially many domains with smaller values of these fields as well, and the properties of 
these domains will not depend on the time t when they are created; see ^U] where this situation 
is discussed. In such models we come to a peculiar conclusion that the main fraction of the 
physical volume of the Universe is in a state with an indefinitely large M p . This might not be a 
real problem, since life of our type simply cannot exist in the parts of the Universe with too large 
(and too small) M p . 

Still it may be important to have global stationarity, see Appendix. The simplest way to 
achieve it would be to impose absorbing or reflecting boundary conditions at sufficiently large 
a, which would preclude the motion of the distribution P p towards large a. Such boundary 
conditions are not unreasonable. Indeed, it is hard to expect that in realistic theories inflation 
will be possible at all indefinitely large values of and a. It may happen, for example, that 
the effective potential V(a) becomes steeper at large a, and inflation (or at least the process 
of self-reproduction of inflationary domains) becomes impossible there. For example, one may 
consider a potential V ~ Xa n e a u . In this theory inflation becomes impossible at a > <Jb = a" 1 . 
As a result, the distribution P p acquires a maximum somewhere near to the Planck boundary 
close to a = 05. Another possibility which leads to a similar effect is that the effective potential 
V(a) decreases at sufficiently large a, for example, V ~ Xa n e~ a a . In this case the distribution 
P p moves to large a until it reaches the maximum of V(a). 
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As we will see from the results of our computer simulations, in both cases the effect of the 
modification of V(a) at large a can be mimicked by the introduction of a boundary at some value 
of the field a = o\>. In our analytical investigation of P p we will assume that 



The way we impose boundary conditions in our numerical investigations will be explained in the 
next section. 



6 Computer Simulations 

The stochastic equations (|37|) are rather complicated partial differential equations, and it is not 
always possible to obtain their solution analytically even in the theories with one scalar field 
[TUj . In the two-field case the situation is even more complicated. Therefore, instead of solving 
these equations directly, we will make a computer simulation of the processes we are trying to 
investigate. 

The main idea of our simulations is the following. We consider N points in the (a, (p) plane. 
Each such point represents the value of the scalar field in a region of size 0{H~ 1 {a ) (p)) (/i-region). 
Our calculations should give us the function P p (a, 0), which is interpreted as the number of h- 
regions with field values a and (p. In our figures this function looks as a two-dimensional surface 
in a three-dimensional space (ex, <p, P p ). 

The values of the fields in each point are initially set to (<jj, (pi). Then we calculate the values 
of the fields in each point independently, since each such point represents an /i-region causally 
disconnected from other h- regions ('no-hair' theorem for de Sitter space). 

Each step of our calculation corresponds to a time change At = uHq 1 , where H = H(ai, (pi), 
and u is some number, u < 1. (The results should not depend on u if it is small enough.) 

The evolution of the fields in each domain consists of several independent parts. First of all, 
each field evolves according to classical equations of motion during inflation. Secondly, each field 



makes quantum jumps by 5a = —yfj^ sinri, 5(p = —yjjj^ sinr2. Here r± and r2 are random 
numbers which are different for each point. 

To make a computer simulation of this branching process, we follow each domain until it 
grows in size two times, and after that we considered it as 8 independent /i-regions. 5 If we would 
continue doing so for a long time, the number of such regions (and our distribution P p (a, (p)) would 
grow exponentially, and it would be extremely difficult to continue the calculation. However, in 
order to obtain a correct probability distribution it is not necessary to follow all domains, since all 
of them evolve absolutely independently. In order to obtain a normalized distribution P p (cr, (p), 
after each step of the calculations we were randomly removing some of the domains, but we were 

5 Note that this does not necessarily correspond to taking steps At = Hq 1 log 2. Indeed, in the domains with 
H(a, </>) 3> Hq the size of the domains grows two times during a time interval much smaller than Hq 1 log 2. 



*(<t,0)L 6 =o. 



(43) 
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doing it is such a way that the probability for any domain to be removed was proportional to 
the distribution P p (a, <fi) at this step of calculation. This allowed us to keep the total number of 
domains fixed and the distribution P p (a, <p) properly normalized without changing at any stage 
of calculation the correct shape of the distribution P p (a, 0). 6 

A special care should be taken about the points near the boundaries. As we already mentioned 
in the previous subsection, there are boundaries of three different types in our problem. 

1) The end of inflation boundary. In the theories with V(a) ~ \a n this boundary is 
given by a ~ ;4=j <fi. When the field a inside a given /i-region becomes smaller than ^= 0, 
inflation in this domain ends, and the value of the field (and of the gravitational constant) in 
this domain practically does not change after that moment. We discard all domains where this 
happens. Then we add new ones in order to preserve correct normalization of our distribution. 
However, each time when we are adding new domains, we distribute them with the probability 
distribution proportional to P p at that time. As we already mentioned, this method allows 
us to keep the distribution P p normalized at all times without distorting its shape. 7 If the 
probability distribution is not stationary, we follow development of our distribution at every 
step of our calculations. However, if the distributions become stationary, one can obtain much 
better statistics by integrating the distribution beginning from the moment when it approaches 
stationary regime. 

2) The Planck boundary. Here we may impose different boundary conditions, depending 
on our assumption concerning the Planck-scale physics. Fortunately, the results which we obtain 
are not terribly sensitive to these assumptions. 

The simplest condition is to discard all points which jump over the Planck boundary, and 
to renormalize the probability distribution P p in the same way as we are doing when the points 
jump over the boundary where inflation ends. 

3) The boundary at large a. In order to obtain a stationary solution we may need to 
have an additional boundary at large a. We will assume that there exists a boundary at some 
sufficiently large value of the field a = For simplicity, we will impose the same condition at 
this boundary as at all other boundaries: we discard all points which jump over this boundary, 
and renormalize the probability distribution P p after such jumps occur. 

An alternative possibility is to consider that the effective potential V(a) becomes very steep 
at large a, and the distribution P p becomes stationary without any need for imposing additional 
boundary conditions at large a. 

6 There is some subtlety here. The volume corresponding to each /i-region is proportional to H~ 3 (<r, 0). Thus, if 
we are interested in the relative fraction of the volume of the Universe, we should show in our figures not the total 
number of /i-regions with given values of the fields, but the total number of such regions multiplied by H~ 3 (<r, <p). 
However, typically the difference between these two distributions is not important, since P p (a,4>) depends on a 
and <p much stronger than H~ 3 (a, </>). 

7 To avoid misunderstandings, we should emphasize that this method cannot not lead to any artificial prolon- 
gation of the stage of inflation. At the stage of self-reproduction of the universe the number of new independent 
domains of the size H^ 1 created due to quantum fluctuations and expansion of the universe is much greater than 
the number of domains disappearing at the boundaries. 
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To plot the distribution we make a two-dimensional histogram of a and and fill the his- 
togram with the points corresponding to each new step of our calculations. After the distribution 
approaches a stationary regime, we instead make a histogram which includes all points starting 
with the step at which the distribution became almost stationary. This does not change the 
shape of the stationary distribution, but effectively increases the number of points involved in 
the calculation, and decreases relative deviation of our 'experimental results' from the probability 
distribution P p . After several hundred more steps we get a rather smooth picture of the stationary 
distribution P p (a, 0). 

In this paper we will present the results of our simulations and of analytic investigation for 
potentials V(a) of several different types. The results of our calculations will be represented as 
a distribution P p inside a box with axes x and y corresponding to the values of the fields a and 
0. The field a grows along the x-axis from a = in the left lower corner. The field grows from 
= when one goes upwards along the y-axis from the same corner. The height z of the surface 
in the box will correspond to the value of P p (a, 0; t). We will not make any attempt to make our 
calculations with realistically small or large values of parameters; our purpose is just to present 
the most important qualitative features of the distribution P p . 



7 Stochastic Processes in Brans— Dicke Theory with a 
Constant Vacuum Energy Density 

In order to get some insight into the complicated behavior of two fluctuating scalar fields, a and 
0, we will temporarily make two steps back to simplify our model. First of all, we will consider 
the theory with the simplest effective potential V(a) = V = const. Also, we will return for 
a moment from the Brans-Dicke theory to the standard Einstein theory. This is equivalent to 
keeping the field = M p ^J~^ fixed. In this case H 2 = Hq = |^ = const. The diffusion 
equation for P p in this theory looks very simple, 

dP p (a,t) _ Hi d 2 P p (a,t) 

~^T~ ~ 8^ da 2 + 3H ° Pp{a > t] ■ (44) 

The solution to the diffusion equation in the comoving frame is a Gaussian with increasing 
dispersion A 2 (t) = J% t. Since the potential is constant, diffusion will always dominate classical 
motion and the Universe will be eternally self-regenerating, with a probability distribution in the 
physical frame P p (a,t) ~ exp(3if ^) exp(— ^A^ft) ) g rows exponentially with time, due to 
the increase in the physical volume of the Universe. 

Let us now consider the same constant potential V(a) = Vq in the Brans-Dicke theory. The 
classical equations of motion read 

(IT 

As before, quantum diffusion of the field a dominates its classical motion, and the Universe will 



(7 = 0, = = const . (45) 
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be in the stage of eternal self-reproduction for the field in the interval 



Let us now analyze the behavior of the probability distribution P p when the field enters 
the self-reproduction range (|4fij). Since H oc 1/0, the distribution will move towards the Planck 
boundary = <p p (|4*U|) along the 0-direction. Very soon the distribution P p reaches the Planck 
boundary and remains concentrated in a very narrow region near it, see Fig. 1. As a result, the 
only possible changes of the distribution P p become related to the Brownian motion of the field a 
along the Planck boundary. The jumps of this field along the Planck boundary are proportional 
to ifp/27r, where H p is the Hubble constant near the Planck boundary, H% = ^ V 1/2 . This 
suggests (and the results of our computer simulations apparently confirm this conjecture) that 
the behavior of P p at large t can be approximately described by equation for one-dimensional 
diffusion along the Planck boundary, similar to eq. 



BP H 3 d 2 P 

The solution to (|47|) is, like in general relativity, a Gaussian with increasing dispersion A 2 (t) = 

H 3 

t and an exponentially growing factor that accounts for the increase in volume at the Planck 
boundary, 

P p (a,t) ~ exp(3H p t) exp L ^^- j . (48) 

It is clear that in this case the distribution gradually becomes flat everywhere along the Planck 
boundary. Therefore the main part of the volume of the Universe will be in a state with indefinitely 
large a. However, as soon as the dispersion A(t) becomes greater than the distance between <To 
and a, the ratio of the volume occupied with a field a to the volume containing field gq becomes 
time-independent. This is an example of the runaway stationarity that we described in the 
Introduction. It is essential that the potential be sufficiently flat in order to have this kind of 
stationarity. For example, such a regime may occur in the theories with V(a) ~ Vq(1 — e _CT ). 

On the other hand, a potential V(a) ~ a 2n will not present this runaway stationarity, since 
any cr-dependence in if oc <r n /0 will make the distribution move forever towards large values of 
u until it is strongly peaked at infinity, unless we impose an extra boundary condition at a = o^. 
However, the behavior of the distribution along the 0-direction will follow the same pattern as 
above, being rapidly concentrated along the Planck boundary (in general a complicated curve, 
depending on the form of the potential V(a)). Its subsequent evolution will be reduced to 
essentially a one-dimensional diffusion along the Planck boundary. 

To illustrate this feature, we performed computer simulation of diffusion for the simplest case 
of the theory with V(a) = const, see Fig. 1. The first few images (Figs, la - lc) correspond 
to the view "from the top" . The horizontal line across the box in Figs, la - lc corresponds to 
the Planck boundary; the distribution P p is concentrated above this line. In the beginning we 
have a delta-functional distribution concentrated near some initial values of a and 0. Then it 
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(c) 



(d) 



Figure 1: Diffusion of the distribution P p along the Planck boundary in the simplest theory with 
V(a) = const, in the time t parametrization. Figs, la, lb, lc show different steps towards stationarity 
at the Planck boundary. Fig. Id shows the same stage of the diffusion as Fig. lc, but in a different 
perspective. 

rapidly moves towards the Planck boundary; it looks like a round spot from the viewpoint we 
have chosen, see Fig. la. After that the distribution widens in the a-direction along the Planck 
boundary, while preserving its width in the ^-direction orthogonal to this boundary, see Figs 
lb, lc. Its shape can be better understood from another viewpoint, see Fig. Id, which shows 
the same distribution as Fig. lc in a different perspective. The most important feature of this 
distribution is that its evolution very soon becomes effectively one- dimensional, being entirely 
concentrated near the Planck boundary. We will take advantage of this feature for the study of 
the runaway solutions in the next Section. 
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8 Runaway Solutions 



Let us study the behavior of the probability distribution P p along the Planck boundary In the 
case of an increasing potential like \a 2n , the larger the value of a in a given domain, the greater 
the increase in physical volume of that domain. Therefore, the probability distribution P p will 
tend to move towards large a. The way it moves will depend on the type of potential. For some, as 
we will see, it is an explosive behavior. The probability distribution gives a statistical description 
of the quantum diffusion process towards large a, but it proves useful to analyze the particular 
behavior of those relatively rare domains in which the field a increases in every quantum jump of 
amplitude H/2n. We can compute the speed at which those domains move towards large values 
of the field a from the equation (where At = H~ l ) 



t _ 5a Hp 
° = At ~ 2tT 




(49) 



where H p is the Hubble parameter along the Planck boundary. For n = 1 there is an exponential 
increase of a in those domains, while for n > 1 we find an explosive solution 



*(t) = °o 1 - -^1^ o n - 1 t . (50) 




We see that for all n > 1 those first domains of the diffusion process reach infinity in finite time. 
Note that the total volume of such domains at that time will be finite, and then they will start 
growing at an infinitely large rate. This behavior is explosive and will correspond to probability 
distributions that are nonstationary and singular at a — > oo. It is extremely difficult to study 
this regime using computer simulations, but with the help of the results obtained in the previous 
section we can get a pretty good understanding of the behavior of P p in such a situation. 

The qualitative analysis and the computer simulations of Section El suggest that the general so- 
lution to the diffusion equation for both fields will factorize naturally into a motion perpendicular 
to the Planck boundary that will reach stationarity very quickly, and a motion along it towards 
large values of a. We will try to study this last motion in the absence of a boundary condition at 
<jfe. We assume that we can ignore the classical motion along the Planck boundary. Indeed, for 
large u and small masses and coupling constants classical motion is almost exactly orthogonal 
to the Planck boundary. This assumption might be violated at large a for some theories with 
very rapidly growing potentials, but in the most interesting case of the theory Ac 4 the classical 
motion is exactly orthogonal to the Planck boundary for all a. On the other hand, for large uj 
and small masses and coupling constants the motion along the Planck boundary almost exactly 
coincides with the motion along the a-axis. In this case an approximate diffusion equation for P p 
analogous to eq. (jUj) can then be written as follows: 

dP p 1 d / r „, /0 d 



H 3 p / 2 —(H^P p )\+3H p P p , (51) 



Ot 8tt 2 da V P da 

where the Hubble parameter along the Planck boundary in the theories with V = i^& 2n is given 
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by 



, 8tt V(a) 8tt 




= Ti# = v --l-l < 52 > 



Equation (|5*T]) can be written as 

dty d 2 ^ 



where 




+ a s p ^ , (53) 



(54) 



4-3n 



Let us first analyze the case V(a) = Acr 4 /4. This corresponds to p = —2. It is clear from (}53|) 
that there are no stationary solutions to this equation in the absence of a boundary condition 
for a. Furthermore, we know from quantum mechanics that potentials of the type — 1/s 2 have 
singular solutions at s = (a = oo). Therefore we expect the distribution P p to be singular at 
a = oo. This was also expected from the analysis of those first domains that reach infinity in 
finite time. 

On the other hand, for the theory V(a) = m 2 a 2 /2, the Schrodinger potential (p)3*j) is of the 
type — s 2 , which has non singular solutions at large a. For small s the solution of eq. (JSlIj) is a 
Gaussian centered at s(t) = (s), with dispersion A(t). The potential then acts asymmetrically 
on it, pulling more at large s and leaving a long exponential tail at small s. The maximum will 
move towards large a, while maintaining a regular solution at infinity. This is expected from the 
previous analysis of the most rapid domains (j4Hj) . For n = 1 we find that the first domains will 
take an infinite time to reach infinity, giving a regular solution at any time. 

However, in both cases we do not have runaway stationarity. Indeed, let us consider equation 
()51|) and assume (for simplicity only) that in the very beginning the function P p was constant. 
Then both for n = 1 and for n = 2 the first term in the r.h.s. of this equation initially is 
positive. Neglecting this term, we obtain P p ~ exp(3if (<r)t), which does not exhibit any runaway 
stationarity. Taking into account the first term in the r.h.s. of equation (J51|) makes the growth 
of P p at large a even faster. This confirms our expectations that in order to obtain runaway 
stationarity one should have an extremely flat effective potential. 



9 Stationary Distributions for Various Theories 

In this section we will study stationary probability distributions for several different potentials 
V(a). Our investigation will mainly rely on the results of our computer simulations, but we 
will try to make analytical investigation whenever possible. As we have argued in the previous 
sections, the simultaneous diffusion of both fields can be approximated by a quick diffusion in 
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the 0-direction towards the Planck boundary and a subsequent diffusion along it until it reaches 
the boundary at er = cr 6 , or until a stationary distribution is established for some other reason. 

1) V(a) = j er 4 . In this case the motion along the Planck boundary is governed by (|53|) with 
p = — 2 and a = 727r/vA and with the boundary conditions ^(0) = \I/(s&) = 0. It is difficult to 
find an exact analytical solution to this simple equation, but one can easily solve it in the WKB 
approximation, 



E — 



a 



Ms) 



where E 
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-1/4 



-1/4 
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Es > a, 



Es < a, 



727TQ-j) 



(55) 

{ r K\f\/Q) 1 ^). This solution has a very sharp maximum close to the 



boundary a = and an exponential decay for small a (large s 2 = 1/er). 
precisely what we observe in the numerical solutions described below. 



This behavior is 



We take the following parameters for our computer simulations: u = 50, A = 0.3. In the 
beginning of the series of calculations we took the points with coordinates (er, <p) close to the 
Planck boundary, but in different initial positions with respect to the boundary ov We have 
found that the duration of the intermediate non-stationary regime depends on the initial values 
of and a. However, typically the stationary distribution P p (a, 0) is established very rapidly, 
after just a few steps At ~ Hq 1 . The resulting stationary distribution is presented in Fig. 2. 
Line OA at this figure corresponds to the end of inflation; there is no inflation for (er, 0) to the 
left of this line. Line OB corresponds to the Planck boundary; V(a) is greater than the Planck 
density M p (0) under this line. The line a = Ob is an additional boundary discussed above. 




Figure 2: Stationary probability distribution P p in the plane (a, 0) for the theory V(a) = fcr 4 , m 
the time t parametrization. The line OA corresponds to the end of inflation, while the line OB to the 
Planck boundary. Inflation occurs between these two lines. The line above OA shows the probability 
distribution V p . 



The curve above the line OA is of the most interest for us. It represents the stationary flow 
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of domains crossing different parts of the boundary OA at the end of inflation. This gives us 
the probability distribution V p (4> e ) that at the end of inflation the field takes some particular 
final value e . The maximum of this curve corresponds to the most probable value of the Planck 
mass M p = \f^-<fi e at the end of inflation. The top of the 'mountain', which shows the stationary 
distribution P p (a,4>), corresponds to the most probable value of M p during inflation. Not unex- 
pectedly, the curve above the line OA looks like a shadow of the mountain P p . Indeed, as we 
have already mentioned, the distribution P p and the distribution V p are directly related to each 
other, see (jHU). Meanwhile, the shape of the distribution P p at small a and large is obviously 
related to its shape at large a and small 0, since in the intermediate regions the points (0, cr) 
follow classical circular trajectories, see e.g. (jHJ)- Consequently, the position of the maximum of 
V p can be approximately obtained by drawing a circle with the center at = o = 0, which goes 
through the crossing point of the Planck boundary and the boundary cr = a b . 

2) V(a) = exp t-. The exponential term is added here in order to show that one can 
avoid introducing additional boundaries at a — a b if the effective potential becomes very steep 
at large a. The result is that the WKB solution (jHSJ), instead of decreasing sharply to a = a b , 
decays exponentially fast. It is just the effect of substituting an infinite barrier by an exponential 
barrier. As we see in the numerical simulations, with a proper choice of the place where the 
effective potential becomes very steep one can reproduce the same result as if there were a 
boundary at a b ~ cr , see Fig. 3. 




Figure 3: Same as in Fig. 2 for the potential V{a) = exp The distribution V v is shown on the 
left side of the box. 

In this figure we show the boundary of the end of inflation by a somewhat wavy line above the 
Planck boundary. At small a the boundary of the inflationary region goes as a straight line from 
the point a = 0, = (compare to the line OA, Fig. 2), but then it becomes curved because 
of the exponential term which precludes inflation at large a. Finally this line crosses the Planck 
boundary. The distribution P p is surrounded by this line and the Planck boundary. 

On the left wall of the box we show the distribution V p , which in the previous picture we 
have shown above the line OA. In other figures we will do the same everywhere when the end of 
inflation boundary is significantly curved. 
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3) V(a) = exp(— s - s ). In this case there is a sharp cut-off of the effective potential at 
o > <Jo, which also leads to the existence of a stationary solution, as if there were a boundary 
near <7o, see Fig. 4. As in Fig. 3, the wavy line corresponds to the end of inflation boundary. 




Figure 4: Same as in Fig. 1 for the potential V(a) 
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-a 2 . In this case the Planck boundary is not a straight line but a parabola 



cr. The eigenvalue equation associated with the probability distribution along the 
Planck boundary is (p)3*j) with p = 2 and a = 144:Tc^/2/m. In the WKB approximation 
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— 1447 ^^ (1 — ( 16^ ) 1 ^ 3 )' This solution has a very sharp maximum close to the 
boundary a = and an exponential decay for small a (small s A = a). This behavior is precisely 
what we observe in the numerical solutions. The distribution P p is very similar to that of the 

A 

1 



theory j a 4 , see Fig. 5. 



5) V(a) 



a' 2 + ^o" 4 log — . Naively, one could expect that the main part of the volume 



TO 



of the Universe in this model should originate as a result of inflation beginning from the points 
on the Planck boundary with the smallest angle 9. Indeed, eq. (JBI) shows that the smaller is the 
initial angle, the greater is the degree of inflation ^7j. However, due to the self- reproduction of 
inflationary domains and the more rapid expansion of domains with greater along the Planck 
boundary, the distribution P p does not stay near the point with the smallest 9, but moves towards 
the largest possible (f> and a, see Fig. 6. As we discussed in the last section, this is a general 
result for any increasing potential, so P p is expected to have a maximum close to the boundary 
a = a b . 

6) V(cr) = -^{jn 2 — Act 2 ) 2 . This is a typical potential used in the theories with spontaneous 
symmetry breaking. It has a minimum at ctq = 77^- in this theory we have two alternative 
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Figure 6: Same as in Fig. 1 for the potential V{a) = ^- a 2 + ^cr 4 log ^. 

regimes. If one begins at o = 0, one may have an inflationary regime at small a ^H], similar to 
the inflationary regime in the new inflationary universe scenario. If, on the other hand, inflation 
begins at large a, then one has an inflationary regime similar to that in the theory ^- a 2 + ^cr 4 . 
The only difference is that in the theory under consideration the field a eventually rolls down not 
to a — 0, but to a = a . The first possibility is illustrated by Fig. 7. On this figure the point 
a = corresponds not to the left corner, as usual, but to the center of the x-axis. 

Note that if one begins with several inflationary domains with different initial conditions, at 
small a and at large a, the domains with large a always win, and the distribution P p very soon 
becomes almost entirely concentrated at a > o"o, see Fig. 8. 
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Figure 7: Same as in Fig. 1 for the potential V(a) = -^{m 2 — Act 2 ) 2 . It describes the stationary 
distribution for the case that inflation begins at a < a Q . 



10 r— parametrization 



The usual Fokker-Planck equation is written in terms of a time parameter t as measured by 
the synchronized clocks of comoving observers. However, in general relativity one can use many 
different time parametrizations. For example, one can measure time by the local growth of the 
scale factor of the Universe and define a new time parameter [2j 

r = In = /* dt' H(a(x, t'), 4>{x, t')) . (57) 

a{x,0) 



This 'time' proves to be rather convenient since in this time, by definition, all parts of the Universe 
expand with the same speed ~ e r , and P p is proportional (though not equal [TO]) to P c . 8 One can 
easily derive the classical equations of motion and quantum diffusion in the new parametrization, 



d(f> _ M 2 (0) dH da M 2 (0) dH 
Or ~ 2irH ~06 ' !h ~ AnH ~da~ 



(5f 



d i ,,\ H 2 2 / , Q \ 9 / 2 \ if 2 2m 2 / 2 v 



dr v / 4tt 2 w v / ' (9r \ / 4tt 2 3iJ 2 

The conditions for the self-reproduction of the inflationary universe are the same as in the 
time t parametrization ([29)1 . Furthermore, we are interested in the diffusion equation for the 

8 Note that dt = F _1 (cj, 4>) dr. Therefore if one is interested in the probability distribution over the invariant 
four-dimensional volume, one should take into account the corresponding sub-exponential corrections to the leading 
e 3r dependence of the 3D-volume on the time r. 
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Figure 8: Same as in Fig. 7 for initial conditions both at large and small a. This series of pictures shows 
that if one begins with equal number of "points" at large and at small a, the volume corresponding to 
large a always dominates in the limit t — > oo. 

scalar fields in the physical frame, where P p (a, <fi; r) = P c (a, 0; r) • e 3r satisfies 

dP L _ d ( M*((j>) dH H d / H 

dr da\4nH da p An da V2tt p // 

, d (M*(4>)dH H d ,H \\ m 

+ d<f>\2irH d<P p+ 4vr^ ( v27r p )) + p ' 

In this case, thanks to the absence of the 3HP p term, there is a candidate for an exact stationary 
solution of eq. (|50|l given by 

P p (a,0;r)cxe 3 ^- 1 (a,0)exp^^^^ . (61) 

Note that this expression is proportional to the square of the Hartle-Hawking wave function of 
the Universe. Unfortunately, this 'stationary solution' does not actually exist in any realistic 
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model of inflation for the reason explained in [10] : The maximum of this distribution coincides 
with the position of the minimum of V(a) where there is no inflation and our stochastic equations 
do not apply. To find a correct solution, we should impose boundary conditions [lOj, 

^(#(a,0)P p (a,0))| cnd = O, 

P P (a,0)|pi anck = O , (62) 
P>,0)| CT6 =O, 

which considerably modifies the shape of the distribution P p . Nevertheless, the naive 'solution' 
floTj) tells us something important about the shape of the distribution P p . We expect this dis- 
tribution to be peaked at the end of inflation boundary, instead of being concentrated near the 

Planck boundary as in the time t-parametrization. The exponent exp ( 3 %y^ 1 remains constant 

all the way along the end of inflation boundary for the theory Xa 4 . Therefore in this theory one 
expects the distribution to move towards small <fi and a, where the prefactor if -1 in ([61]) is max- 
imal. In other theories this argument does not apply, and in general one may obtain a stationary 
distribution with parameters depending not only on the boundary of the end of inflation but also 
on the boundary at a/,. All these features are observed in the numerical simulations. 



11 Computer Simulations of diffusion in time r 

Computer simulations in the time r are similar to the ones in the time t, but there are several 
important differences. 

First of all, there is no need to make any splits of the domains. The reason is that all domains 
in the time r expand with the same speed, by definition of this time as a logarithm of expansion. 
Each step of our calculation now corresponds to a time change At = u, where u is some small 
number, u < 1. 

As before, evolution of the fields in each domain consists of several independent parts. Each 
field evolves according to classical equations of motion, (|58|). and in addition to it each field makes 
quantum jumps by 8a = sinr 1; 5(p = sinr 2 . Here r 1 and r 2 are random numbers which 
are different for each point. 

As in Section 03 in addition to the distribution P p ((p, a, r) we computed also the volume of all 
domains where inflation ended within each new time interval At. This gives the fraction of the 
volume of the Universe where inflation ends at a given time r within a given interval of values of 
the field 0. We call this distribution P p (0 e ,r). In the theories with V(a) = ^cr 2n 

V p ((j)e, t) ~ (j) e ■ P p ((fre, cr e , t) ■ (63) 

It is very instructive to compare the results of computer simulations in time t and in time r 
for various effective potentials V(a). The general tendency we observe is that the distribution 
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P p (0, a, t) is more closely concentrated not near the Planck boundary or near the boundary at 
large a (as was the case for the distribution P p ((f>,a,t)), but near the boundary corresponding 
to the end of inflation. The reason for this difference is very simple. In the time r there is no 
additional enhancement of the volume filled by the fields corresponding to large values of the 
Hubble constant. Now let us consider several particular examples. 

1) V(a) = const. This corresponds to the simple model we considered in Section 6. In the 
time t the distribution rapidly moved towards the Planck boundary and then diffused along it, 
Fig. 1. Evolution of P p in the time r is quite different, for the reason discussed above. The 
distribution moves away from the Planck boundary, see Fig. 9. This explains many features of 
the distribution P p in more complicated theories to be considered below. 




Figure 9: In the theory with V{a) = const the distribution P p (4>, c;r) moves away from the Planck 
boundary. Compare to Fig. Id. 

2) V(a) = j a 4 . In this theory we do not have any typical mass scale which would correspond 
to a maximum of the probability distribution. Therefore in time t the distribution was moving 
towards large a and <ft, until it was stabilized either by a boundary at large a or by the change 
of the potential at large a (e.g. V(a) = \a 4 ■ exp^r, see next item). In the present case the 
distribution moves towards smaller and smaller a and 0, and there is no stationary regime unless 
the potential V(a) at small a becomes, for example, quadratic in a, see below. Fig. 10 shows 
how the distribution moves towards the corner a = <fi = 0. This Figure should be compared to 
Fig. 2, where the corresponding distribution is shown in time t. 

As in the Section |5J the curve above the line corresponding to the end of inflation shows the 
probability distribution V p (<p e , r) that at the end of inflation the field takes some particular 
final value e . Note that in the present case this distribution depends on r and moves towards 

0e = O 

3) Via) = -jcr 4 ■ exp q -^. No qualitative difference appears here as compared with the theory 
| a 4 since the exponential term does not modify the potential at small a. 

4) V{a) = \<j 2 . In this case the distribution P p is stationary, see Fig. 11. As one can see, 
it is concentrated near the boundary of the end of inflation. This Figure should be compared to 
Fig. 5. 
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Figure 10: Probability distribution in the plane (cr,(f>) for the theory V(a) = fo" 4 , in the time r 
parametrization. Compare to Fig. 2. 




Figure 11: Same as in Fig. 10 for the theory V(a) = ^a 2 . 

5) V{a) = ^a 2 + | <x 4 log ^. In this model, and in a simpler model with V{a) = ^ a 2 + \a\ 
the quadratic term stabilizes the distribution P p , see Fig. 12. The resulting distribution is 
stationary even in the absence of the boundary at large a; compare with Fig. 6. 

6) V(cr) = -^{m 2 — Act 2 ) 2 . Here we have two alternative regimes. If one begins at large a, 
then one has an inflationary regime similar to that in the theory ^ a 2 + \& A . The corresponding 
distribution will be very similar to that shown in Fig. 12. One may have an inflationary regime 
at small a |lfij . similar to the inflationary regime in the new inflationary universe scenario. The 
corresponding distribution is shown in Fig. 13; compare with Fig. 7. 

Note that the stationary distributions we have obtained look different from those for the same 
theories in time t. This means that the normalized distribution P p in the large time limit does 
not depend on time, but it does depend on the choice between different 'times' (t, r, etc.) A 
similar conclusion was earlier reached for other models studied in JTU] • The reason of this strange 
behavior can be understood by using the following simple analogy. 
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Figure 12: Same as in Fig. 10 for the potential V(a) = \ a 2 + f cr 4 log ^. 




Figure 13: Same as in Fig. 10 for the potential V(a) = -^{m 2 — Xa 2 ) 2 . It describes the stationary 
distribution for the case that inflation begins at a < a Q . 

Let us consider a two-dimensional plane (x, y) and a cone formed by the lines y = x and 
y = — x going from the point (0, 0) towards positive y, see Fig. 14. Let us paint light grey the 
area inside this cone to the right of the ?/-axis, and paint gray the area inside the cone to the left 
of the y-axis. Now we will consider y as a time direction, cross the cone by the lines of fixed y at 
a distance dy from each other and compare the white area and the gray area in the interval from 
y to y + dy. Obviously, for all dy the ratio of the white area to the gray area will be equal to 1, 
independently of the time y (stationarity). 

Now one may choose another 'time' direction y' by rotating the y-axis, and slice the cone by 
the lines of constant y' and y' + dy'. In this case the ratio of the white area to the gray area in the 
interval from y' to y' + dy' also will not depend on the time y' (stationarity), but this ratio will 
not be equal to 1 anymore. This strange effect is possible due to the fact that the total area of 
the gray part of the cone, as well as of the white one, is infinite. As usual, when the integrals are 
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Figure 14: A cone in a two-dimensional space formed by the lines y = x and y = —x. Two different 
parametrizations of "time" y give different ratios for the gray over the white area, although they are 
both stationary. 

divergent, their ratio depends on the way one takes them. A similar effect appears in our case as 
well. The total volume of all inflationary domains in a self-reproducing universe is infinite in the 
limit t — > oo (or r — > oo). Therefore the relative fraction of the volume of all domains with any 
particular properties may depend on the way in which we are sorting out these domains. This is 
the main reason of the difference between P p (a,<f);t) and P p (cr, 0;t). This difference exists even 
if each distribution is stationary. 



12 Discussion 



Let us try to summarize our results and discuss their possible implications. First of all, we have 
confirmed that the regime of self-reproduction is possible not only in the ordinary inflationary 
theory, but in the Brans-Dicke inflation as well. However, the Brans-Dicke inflation has some 
new interesting features. In this theory the upper (Planck) boundary for the energy density of 
a classical space-time is not a point V(a) = as in the Einstein theory, but a line V(a) = 

^r<t> A - Similarly, the end of inflation boundary in the simplest models of chaotic inflation with 
V(a) ~ a n is not a point o ~ M p but a line a ~ <p/^/uj. In the ordinary inflationary models the 
probability distribution P p to find the inflaton field a at a given time in a given volume typically 
is concentrated either near the Planck boundary or near the boundary where inflation ends 1CJ. 
In the Brans-Dicke theory with V(a) ~ o 11 or with V(a) ~ e acr this distribution also approaches 
one of these two boundaries, but after that it may continue moving, sliding along the boundaries. 

As a result, the probability distribution P p approaches stationary regime only if there exist 
some additional reasons which preclude this sliding. This may happen, for example, if the effective 
potential becomes very steep (or if it decreases) at large a. We have studied this possibility both 
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by making computer simulations in the theories with potentials which become rapidly increasing 
(or decreasing) at large a, and by introducing a phenomenological boundary at a = er 6 . We 
obtained stationary probability distributions for a wide class of theories with two different time 
parametrizations. 

If nothing precludes sliding of the probability distribution along the boundaries, one typically 
obtains a non- stationary distribution. However, the local stationarity still exists: The properties 
of domains with given values of scalar fields <p an d cr do not depend on the time when these 
domains were formed. 

Whether the probability distribution P p is stationary or not, in the process of its evolution it 
probes all values of the fields and a for which inflation and self-reproduction of the Universe 
can take place. As a result, after inflation the Universe becomes divided into many exponentially 
large domains with different values of the effective Planck mass M p (</>). 

It would be natural to assume that the probability to live in a typical part of our Universe is 
proportional (though not equal, see below) to V p (4> e ). In particular, if our calculations would give 
us a delta-functional distribution V p (<p e ), we would have a definite prediction for M p and thus for 
G. Our results show, however, that the distributions which we obtain for large values of masses 
and coupling constants are rather smooth. If one takes very small masses and coupling constants, 
which is necessary to obtain small density perturbations Sp/p, the distributions become very 
sharply peaked indeed, but still they are not delta-functional. 

The probability to live in a given part of the Universe depends not only on its volume (which 
is proportional to the distributions P p and V p if they are sufficiently narrow) but also on the 
conditions inside this volume. These conditions depend very strongly on the value of M p . For 
example, it is well known that a decrease of the Planck mass M p by less than an order of magnitude 
from its present value in our part of the Universe would make the lifetime of the Sun so small 
that no biological molecules would appear on the Earth. An even bigger decrease of M p would 
lead to an extremely efficient nucleosynthesis and to the absence of hydrogen in the Universe ■ 
An increase of M p would slow the expansion of the Universe. In such a Universe the departure 
from thermal equilibrium during the process of baryogenesis would be small, this process would 
be inefficient and the Universe now would be practically empty. On the other hand, a decrease 
of M p decreases the reheating temperature after inflation, which may also cause the absence of 
baryons. Finally, in the simplest model of inflation with V(a) = ^cr 2 , density perturbations 
produced during inflation are inversely proportional to M p , see eq. (|28j) . Therefore a change of 
M p would lead to a profound modification of the properties of galaxies. 

This suggests that the knowledge of the distribution V p being complemented with anthropic 
considerations may help us determine the most probable value of the gravitational constant in 
the domains of the Universe where life of our type is possible. This is a very exciting possibility, 
resembling the 'big fix' paradigm of the baby universe theory j57j, but, just like the baby universe 
theory, it involves many speculations. One of the problems of such an approach is the dependence 
of V p on the choice of time parametrization. We will briefly discuss this issue in the Appendix. 
Whether this most ambitious part of our program will be successful or not, it is certainly true 
that the theory of a self-reproducing inflationary Brans-Dicke universe offers us many interesting 
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and unexpected possibilities. 

For example, in the standard inflationary cosmology the Planck mass was fixed, and in order to 
obtain a desirable amplitude of density perturbations ^ ~ 5 • 1CT 5 in the theory ^-cr 2 one should 
introduce into the theory a new mass scale, m ~ 10 13 GeV, which is six orders of magnitude 
smaller that the Planck mass M p ~ 10 19 GeV. This cannot be considered as a fine tuning; after 
all, the electron mass is 22 orders of magnitude smaller than the Planck mass. Still it is very 
interesting to see how the same issue looks in the context of the inflationary Brans-Dicke theory. 

In this theory m is fixed but the Planck mass is not. It takes all its possible values in 
different parts of inflationary universe. Correspondingly, the amplitude of density perturbations 
in different parts of the Universe takes all possible values from to 0(1). If by calculating V p 
and using anthropic considerations we are able to explain why we live in a part of the Universe 
with & ~ 5 ■ 1(T 5 (see Appendix), then we will simultaneously explain why is so small and 
why M p is so large in our part of the Universe. But there is a chance that we live in a part of the 
universe with ^ ~ 5- 10~ 5 without any special reason, just as some people live in exotic countries 
without even knowing that they are exotic. Nevertheless, even in this case we will get something 
interesting. We should simply look at equation (|2*%|) in a different way, writing it as follows: 

M p (0) ~ 50 m • (^j . (64) 

This equation tells us that if we live in a part of the Universe with ^ ~ 5 • 10~ 5 , then in this 
part of the Universe the effective Planck mass automatically happens to be one million times 
larger than m, which is the only mass scale we have in our simple theory. In other words, we 
do not need to introduce into the theory two different mass scales different from each other by 
a factor of 10 6 . It is enough to have one mass scale; the rest of the job will be accomplished 
by quantum fluctuations in the inflationary universe. In this scenario both smallness of the 
density perturbations and greatness of the Planck mass appear as two sides of the same purely 
environmental effect. 
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Appendix. Towards determination of the most probable 
value of the gravitational constant G(4>) 



One of the most interesting problems of elementary particle physics is to understand why the 
gravitational constant is so small or equivalently the Planck mass M p is so large. In the context 
of our model this problem is formulated in a very unusual way. Our Universe consists of many 
exponentially large domains with different values of G — M~ 2 (</)). Therefore instead of finding 
a unique value of M p for the whole Universe, one should study the distribution of all possible 
values of M p in our Universe. This might give us a possibility to understand why the gravitational 
constant is so small in the part of the Universe where we live. 

First of all, one can make a natural assumption that the number of observers asking questions 
about the parts of the Universe with given values of fields and a is proportional to the volume 
of these parts. This suggests that the answer to our questions may be related to the investigation 
of the distributions P p and V p . However, it is not clear whether it is possible to justify this 
suggestion, especially if one recalls that these distributions depend on the choice of the time 
parametrization. 

Note, that these probability distributions have a very well determined operational meaning 
when one uses them to predict the distribution of a scalar field in the Universe at a specific 
hypersurface t = const (or r = const) under given initial conditions at t = (or at r = 0). 
Now the main question is whether it is possible to use the results of our calculation of P p and 
V p to get some information about the most probable value of the gravitational constant in those 
parts of the Universe where we can make observations and ask questions about the gravitational 
constant. 

One may argue that as far as we have a nonvanishing probability to live in the domains with a 
given value of M p (</>), and as far as the total volume of all such domains integrated over all times 
is infinite (the integral exponentially diverges at t — > oo), there is no reason to compare these 
infinities and study details of behavior of P p . We live, which means that we have picked up one 
of these domains, but all attempts to go any further than that and to explain 'our choice' would 
not make any sense. This would be similar to attempts of a man from Switzerland to understand 
why he was born there rather than in China where the total number of people is much greater. 

According to this point of view, one should use our results in a very limited way. One should 
find those domains where the distribution V p and the probability of existence of life of our type do 
not vanish. (The answer to this question does not depend on the choice of time parametrization.) 
One should then find out what is the value of the effective Planck mass in each of these domains, 
and to relate it to other properties of space and matter in these domains. 

In the Discussion we took this most conservative attitude towards our results. Even in this 
case inflationary cosmology goes far beyond the standard Big Bang paradigm, which assumes 
that the gravitational constant should be the same in all parts of our Universe. 

This conservative attitude is quite legitimate. It is quite possible indeed that one should not 
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ask why he or she was born in this or that country and in this or that part of the Universe. 
However, this idea is suspiciously similar to the old belief that it does not make any sense to 
question initial conditions in the Universe: The Universe is big and flat for the reason that it was 
born big and flat; it contains more baryons than antibaryons for the reason that it was created 
that way... After the invention of the theory of baryogenesis and of inflationary cosmology this 
way of avoiding complicated problems does not look particularly attractive. 

On the other hand, we are not well prepared to go beyond this point. In order to understand 
why do we live in a part of the Universe with a given value of the gravitational constant, one 
should learn first what is life, how it appears, is it correct that the probability for life to appear 
is proportional to the space available, etc. Before we do it, there will be no guarantee that we are 
on the right track. However, instead of waiting until the theory of everything is constructed, we 
may use inflationary cosmology as a good playground where we can test various hypotheses. This 
might help us to learn how to formulate correct questions (and, if we are lucky, to get correct 
answers) in the context of the new cosmological paradigm. 

To illustrate some ambiguities involved in this approach, let us consider a problem formulated 
some time ago by Holger Nielsen j2H]. Assume that we live in a peak of probability to be 
born at some particular time t. This hypothesis at the first glance looks quite reasonable and 
innocuous. In any case, it does not look obviously wrong if we are looking to the past. Indeed, 
the total population of the Earth now is much larger than it was before, and it continues growing 
exponentially. Most of the people who have ever lived on the Earth were born in the 20th century. 
This gives us a total of less than 20 billion people. If we assume that a typical person should 
be born near the maximum of the probability distribution, and if we assume that we are typical, 
then we would not expect much more that 20 billion people to be born from now on. But this is 
possible only if very soon, within the next few decades, the population of Earth will start rapidly 
decreasing [2T?j . 

We definitely want this doomsday prediction to be wrong, but what could be wrong about 
it? The point is that this prediction is a consequence of the assumption that the total number 
of people to be born is finite, and we live near the maximum of the distribution. However, if the 
total number of the people to be born is infinite, then we live at some time t not for the reason 
that this time is near the maximum of the distribution, but for the only reason that we must pick 
up some finite time t rather than the time t — oo. For example, in our scenario the total volume 
of the Universe grows exponentially at all times. Consequently, the total number of domains of 
our type, the total number of planets of the type of the Earth and the total number of people 
populating these planets also grow exponentially [TU] . Thus, in this scenario the total population 
of the Universe does not have any maximum and any falloff at large t. 

This example shows that one should be extremely careful with probabilistic arguments to 
avoid many hidden ambiguities. For this reason we removed this discussion from the main body 
of the paper, to make sure that our speculations do not mix up with reliable results. 

One of the obvious problems with the distributions P p and V p is the choice of the time 
parametrization. Indeed, as we have seen, position of the peaks of these distributions does 
depend on this choice. Of course, one may argue that life as we know it uses ordinary time t 
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related to periodic processes, rather than the time r which measures the logarithm of the distance 
between galaxies. Does it suggest that the proper distribution to study is P p (a, <fi; t)l Let us make 
this assumption for a moment and see what we will learn. 

First of all, the peaks of the distributions P p for all theories we studied correspond to the state 
with a maximal possible value of the Hubble constant compatible with the self-reproduction of 
inflationary domains. For all theories V(a) ~ a n the peak of P p appears near the upper boundary 
(Jb, and the peak of V v appears near ~ cr fe y^, which gives M p ~ 2<jfty^. The boundary cr& may 
correspond either to the place where V(a) changes its shape and becomes very steep, or to the 
place where it becomes decreasing. 

Now one should multiply the volume of the domains with given M p by the probability that 
life of our type can exist in these domains. The results look rather strange. Anthropic arguments 
seem to allow variations of M p at least by one order of magnitude. Meanwhile typical distributions 
V p are extremely sharp if masses and coupling constants are realistically small. Multiplication 
of these distributions typically gives a distribution which is sharply concentrated near one of the 
boundaries of the anthropically allowed region. This is very similar to the results obtained by 
Rubakov and Shaposhnikov when they discussed anthropic considerations in the context of the 
baby universe theory [3*0] . 

The situation becomes even more complicated when one considers models where the distribu- 
tions P p and V p are non-stationary. In such models the probability distribution V p permanently 
moves towards infinitely large <fi. This growth eventually outweighs all anthropic bounds and 
pushes the distribution towards the region where life of our type could exist only as an extremely 
unstable phenomenon. This would be clearly incompatible with the results of anthropic consid- 
erations, which suggest that we live not far away from the center of the anthropically allowed 
region of possible values of M p . 

One may try to interpret this conclusion as an argument against Brans-Dicke cosmology, 9 
or at least against those versions which lead to the runaway solutions for P p and V p . Another 
possibility is to use the distribution in the time r, which is typically much less sharp. Note 
that this time parametrization also has certain advantages, see [22]. Still another possibility is 
that the distributions P p and V p , being very useful for the description of the global structure of 
inflationary universe at a given time, cannot be used for the calculation of the probability of life 
appearing in the part of the Universe of our type. Indeed, as we know from j2H], the notion of 
time (either t or r) makes sense in the context of quantum cosmology only after the appearance 
of observers. 

It is excitingly interesting to participate in the investigation of this problem. However, at 
the moment we do not even know whether physics provides a wide enough framework to study 
appearance of life, or some additional ingredients are needed pQ. This is an extremely complicated 
and speculative issue. We hope to return to its discussion in a separate publication. In the 
meantime we decided to restrict ourselves to the conservative approach outlined in the Discussion. 



9 There exist some arguments in the context of the baby-universe theory that the Brans-Dicke theory with an 
account taken of the wormhole effects is reduced to the standard Einstein theory |3J . 
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